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^^ I Abstract. We present a two- fluid magnetohydrodynamics (MHD) model of quasi- 

^^ ' stationary, two-dimensional magnetic reconnection in an incompressible plasma 

composed of electrons and ions. We find two distinct regimes of slow and fast 
P^ I reconnection. The presence of these two regimes can provide a possible explanation for 

C/j ■ the initial slow build up and subsequent rapid release of magnetic energy frequently 

(~| , observed in cosmic and laboratory plasmas. 
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1. Introduction 

Magnetic reconnection is the physical process by means of which magnetic field fines 
join one anotfier and rearrange tfieir topology. Magnetic reconnection is believed to be 
tfie mecfianism by wfiicfi magnetic energy is converted into kinetic and tfiermal energy 
in tfie solar atmospfiere, tfie Eartfi's magnetospfiere, and in laboratory plasmas [H El 
[31 m m [HI [7] . Many reconnection related pfiysical pfienomena observed in cosmic and 
laboratory plasmas exfiibit a two-stage befiavior. During tfie first stage, magnetic energy 
is slowly built up and stored in tfie system witfi relatively little reconnection occurring. 
Tfie second stage is cfiaracterized by a sudden and rapid release of tfie accumulated 
magnetic energy due to a fast reconnection process. For example, a solar flare is powered 
by a sudden (on timescale ranging from minutes to tens of minutes) release of magnetic 
energy stored in tfie upper solar atmospfiere [1]. Because tfie value of tfie Spitzer 
electrical resistivity is very low in fiot plasmas, magnetic energy release rates predicted 
by a simple single-fluid MHD description of magnetic reconnection are mucfi slower 
tfian tfie rates observed during fast reconnection events in astropfiysical and laboratory 
plasmas [H El HI |5l |6l [Zj. One of tfie most promising solutions of tfiis discrepancy is 
tfie two-fluid MHD tfieoretical approacfi to magnetic reconnection [H IH El El [71 and 
references tfierein]. Recently a model of two fluid reconnection in a electron-proton 
plasma was presented in [8]. In tfiis paper, we consider a more general case of two-fluid 
reconnection in electron-ion and electron-positron plasmas, and we present derivations 
in detail. In tfie discussion section, we also argue tfiat tfie slow and fast reconnection 
regimes predicted by our model, can provide a possible explanation for tfie observed 
two-stage reconnection befiavior. 

2. Two-fluid MHD equations 

In tfiis study, we use pfiysical units in wfiicfi tfie speed of ligfit c and four times tt 
are replaced by unity, c = 1 and 47r = 1. To rewrite our equations in tfie Gaussian 
centimeter-gram-second (CGS) units, one needs to make tfie following substitutions: 
magnetic field B — )■ B/vivr, electric field E — t- cE/yivr, electric current j — t- yivr j/c, 
electrical resistivity r] — )■ 'qc^/An, tfie proton electric cfiarge e — )■ A/ivre/c. 

We consider an incompressible two-component plasma, composed of electrons and 
ions. We assume tfie plasma is non-relativistic and, tfierefore, quasi-neutral. Tfie ions 
are assumed to fiave mass rrii and electric cfiarge Ze, wfiile tfie electrons fiave mass me 
and cfiarge — e. Because of incompressibility, tfie electron and ion number densities are 
constant, 

Ue = n = const, rij = Z^^n = const, (1) 

wfiere tfie last formula follows from tfie plasma quasi- neutrality condition Zerii = erie. 
Tfie plasma density p, tfie electric current j and tfie plasma (center-of-mass) velocity V 
are 

p = miUi + meUe = n^Z^^nii + nie) = const, (2) 
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j = ZeUiU^ — eUeU'^ = ne(u* — u'^), (3) 

V = [niiniU.'' + rrieneVL^) / p = n{Z~^miVi + me\y^)/ p. (4) 

Here u*^ and u* are the mean electron and ion velocities, which can be found from the 
above equations, 

u^ = V - (mi/Zep)j, u' = V + {mjep)]. (5) 

The equations of motion for the electrons and ions are [9l ITU] 

Uerue [dt\f + (u"V)u"] = - VPe - riee(E + u^ X B) - K, (6) 

Uinii [dtu' + (u^V)u^] = - VPi + niZe(E + u^ x B) + K, (7) 

where Pe and Pi are the electron and ion pressure tensors, and K is the resistive frictional 
force due to electron-ion collisions. Force K can be approximated as [9l [10] 

K = n'^e^r]{u'' - u*) = -ner]i, (8) 

where r] is the electrical resistivity, and we use equation ([3]). For simplicity, we assume 
isotropic resistivity, and we also neglect ion-ion and electron-electron collisions and the 
corresponding viscous forces. Substituting equations ([I]), (I5l) and ([H]) into equations ([6]) 
and ([7j), we obtain 

nme [dtV + (VV)V] 

-{nm,m,/Zep) [dt] + (VV)j + (jV)V - (m,/Zep) (jV)j] 
= —VPe — neE — neV x B + (min/Zp)] x B + nerj}, (9) 

Z-^nmi [dtV + (VV) V] 

+ {nmemi/Zep) [dti + (VV)j + (jV)V + {mJep){jV)}] 
= —VPi + neE + neV x B + {711^12/ p)} x B — ner]}. (10) 

We sum equations iQ and fITOl) together and obtain the plasma momentum equation 

p[dtV+{VV)V] + Km,/ZeV)(jV)j = -VP + j x B, (11) 

where P = Pe + Pj is the total pressure. Next we subtract equation (ITOl) multiplied by 
ZiTLe/mi from equation iQ and obtain the generalized Ohm's law 

E = ?7J — VxB + {mi/Zep){l — Zme/rrii)} x B 
- im,/Zep)[VPe - {ZmJmi)VPi] 
+ {m,m.jZe^p) [da + (VV)j + (jV)V 

-(m,/Zep)(l -Zme/m,) (jV)j]. (12) 

It is convenient to introduce the ion and electron inertial lengths 
d, = (m,MZ2e2)i/2 = (m,/Zne^y/\ 
4 = (me/nee2)i/2 = (mjne'y/' < d„ ^ ' 

and constants 

ujI = {1 + Zmjmi)-^ = (1 + dl/dfr\ 



uo 
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1 - Znie/nii = 1 - dl/df > 0. 



(14) 
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Here we consider a physically relevant case of Znie < nii, so that d^, < di,0 < u'^ < 1 and 
1/2 < w^ < 1. Note that w^ ^ cj^ ^ 1 in the case of electron-ion plasma {Zirie <^ rui), 
and u"^ = 1/2 and w^ = in the case of electron-positron plasma {Z = 1 and mi = me). 
Using definitions (TT3|) and ([HD, we obtain for the plasma density ([2]) expression 

p = m^iTi/Zuj']^ = n'^e'^d'^/u'^, (15) 

and we rewrite the plasma momentum equation (TTTl) and Ohm's law ( iT2l) as 

p[dtV+{VV)V]+uldl{3V)i = -VP + j X B, (16) 

E = r7J-VxB + (ulujl/ne) j x B 

-iul/ne)[\/Pe~idl/d^)VP;\ 

+ uldl [dt] + (VV)j + (jV)V - (a;^a;!/ne)(jV)j] . (17) 

It is noteworthy that the electron inertia terms, proportional to d1, enter both Ohm's 
law and the momentum equation. Although these terms are important for fast two- 
fiuid reconnection (as we shall see below), they have been frequently neglected in the 
momentum equation in the past |j. In addition, we note that V ■ B = 0, and also 
V ■ V = and V ■ j = for incompressible and non-relativistic plasmas. 

For convenience of the presentation, below we will refer to the plasma as being 
electron-ion, even though, unless otherwise stated, our derivations in the next two 
sections are valid for reconnection in an electron-positron plasma as well. 

3. Reconnection layer 

We consider two-fiuid magnetic reconnection in the classical two-dimensional Sweet- 
Parker- Pet schek geometry, which is shown in figure [H The reconnection layer is in the 
x-y plane with the x- and y-axes perpendicular to and along the reconnection layer 
respectively. The z derivatives of all physical quantities are zero. 

The approximate thickness of the reconnection current layer is 25, which is defined 
in terms of the out-of-plane current (j^) profile across the layer |§|. The approximate 
length of the out-of-plane current (j^) profile along the layer is defined as 2L. Outside 
the reconnection current layer the electric currents are weak, the electron inertia is 
negligible. Ohm's law ( [T71) reduces to E = —V x B -|- j x B/ne = —vl^ x B (in the 
case of electron- ion plasma, w^ ~ wi ~ 1), and, therefore, the magnetic field lines are 
frozen into the electron fiuid. Thus, 25 and 2L are also approximately the thickness and 

X For particle species s S {e, i] we use the standard definition of the pressure tensor as the 
density times the second moment of the particles velocity fluctuations relative to the mean velocity, 
Ps = nsms{{v'^ — \x^){v'^ — u")) , where u"* — {v'^) [9^. Instead, one could use velocity fluctuations relative 
to the plasma center-of-mass velocity ^ and define pressure as P, = nsmsdv" — V)(u* — V)) [10]. In 
this case, the total pressure tensor would he P = Pe + Pi = P + w^dg j j, and, therefore, the electron 
inertia term w^(ig(jV)j in the momentum equation (|16|) would become absorbed into the pressure term 
VP. However, note that pressure P is strongly anisotropic. 

§ Thickness 5 can be formally defined by fitting the Harris sheet profile {Bext/S)cosh^'^{x/S) to the 
current profile jzix, y — 0). 
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Figure 1. Two-dimensional reconnection layer. 



the length of the electron layer, where electron inertia is important and the electrons 
are decoupled from the field lines. The ion layer, where the ions are decoupled from 
the field lines, is assumed to have thickness 2 A and length 2Lext, which can be much 
larger than 26 and 2L respectively. The values of the reconnecting field in the upstream 
regions outside the electron layer (at x ~ 5) and outside the ion layer (at x ~ A) 
are about the same. By ^ B^xt up to a factor of order unity. This result follows 
directly from the definition of 26, and from the z-component of the Ampere's law, 
By{x,y = 0) = j^ jz{x',y = 0)dx'. The out-of-plane field B^ is assumed to have a 
quadrupole structure (see figure [1]) O [6l [7] [|jl. 

The reconnection layer is assumed to have a point symmetry with respect to its 
geometric center O (see figure [1]) and reflection symmetries with respect to the x- 
and y-axes. Thus, the x-, y- and 2;-components of V, B and j have the following 
symmetries: Vxi±x,Ty) = ±Vx{x,y), Vy{±x,Ty) = TVy{x,y), V;(±x, =Fy) = Vz{x,y), 
Bxi±x,Ty) = TBxix,y), By{±x,Ty) = ±By{x,y), B^(±x, =Fy) = -5^(x,?/), 
jxi±x,Ty) = ±jxix,y), jy{±x,Ty) = Tjyix,y) and j^(±x, =Fy) = jzix,y). The 
derivations below extensively exploit these symmetries and are similar to the derivations 
inPIIIlig. 

We make the following assumptions for the reconnection process. First, resistivity 
rj is assumed to be constant and very small, so that the characteristic Lundquist number 
S is very large, 

S = VaL^xiN > 1, Va = B^xtl^p. (18) 

Here Va is the Alfven velocity. Second, the reconnection process is assumed to be quasi- 
stationary (or stationary), so that we can neglect time derivatives in the equations 
above and in the derivations below. This assumption is satisfied if there are no plasma 
instabilities in the reconnection layer, and the reconnection rate is slow sub-Alfvenic, 
Eiz "^ VABext- Third, we assume that the reconnection layer is thin, 6 <^ L and 
A ^ Lext, which is an assumption related to the previous one. Fourth, we assume that 
the electron and ion pressure tensors Pe and Pi are isotropic, therefore, the pressure 
terms in equations (fTTl) and (fT6|) are assumed to be scalars. 



II Below we shall see that Bz has quadrupole structure only in the case of electron-ion plasma, but 
not in the case of electron-positron plasma. 
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4. Two-fluid reconnection equations 

We use Ampere's law and neglect the displacement current in a non-relativistic plasma 
to find the components of the electric current 

j^ = dyB^, jy = -d^B^, j^ = d^By - dyB^. (19) 

The ^-component of the current at the central point O (see figure H]) is 

jo = {jz)o = {dxBy - dyB.^)o ^ {d^By)o ^ Bextl^, (20) 

where we use the estimates {dyB^)o ^ {dxBy)o and {dxBy)o ~ B^xt/^ at the point O. 
The last estimate follows directly from the definition of 6 as being the half-thickness of 
the out-of-plane current profile across the reconnection layer. 

In the case of a quasi- stationary two-dimensional reconnection, we neglect time 
derivatives, and Faraday's law V x E = — t^jB for the x- and y-components of the 
magnetic field results in equations dyEz = —dtB^ = and dxEz = dtBy = 0. Therefore, 
Ez is constant in space, and from the z-component of the generalized Ohm's law (ITTll 
we obtain 

Ez = VJz - VxBy + VyBx + {uj'\uj'i/ne){jxBy - jyB^) 

+ Uldl [Vxdxjz + Vydyjz + JxOxVz + jj/^y V^ 

-{uj'Xuj'i/ne){jxdxjz+ jydyjz)] = constant. (21) 

The reconnection rate is determined by the value of Ez at the central point O, that is 

Ez = Wo- (22) 

We see that the electric field is balanced only by the resistive term rjjo at the central 
point O; this is because we assume isotropic pressure tensors in this study. To estimate 
joi in what follows we neglect time derivatives for a quasi-stationary reconnection and 
we use the symmetries of the reconnection layer. 

The z-component of the momentum equation ( !T6|) is 

p{VxdxVz + VydyVz) + (^Idlijxdxjz + jydyjz) = jxBy - jyBx- 

Taking the second derivatives of this equation with respect to x and y at the point O, 
we obtain 

p{dxVx)oidxxVz)o + ujldl{dxjx)oidxxJz)o = {dx3x)o{dxBy)o, 

P{dyVy)o{dyyVz)o + (^id^dy j y) o^Oyy J z) O = - {Oy ] y) o^Oy B x) O ■ 



Therefore, 

{,dxxVz)o = -idxyBz)o[{dxBy)o - ujldl{dxxjz)o\/ pidyVy)o, 
{dyyVz)o = {dxyBz)o[{dyBx)o + ujldl{dyyjz)o]/ p{dyVy)o, 

where we use equations ( fT9l) and the plasma incompressibihty relation dxVx = —dyVy 



(23) 
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Next, we calculate the second derivatives of equation (ET|) with respect to x and y 
at the central point O and obtain 

= v{d..Jz)o - WxV.)o - {culcu^_/ne){d.,j.,)o]{d,By)o 
+ 2ujldl[{dMo{dxxJz)o + {dxJx)o{dxxVz)o 

- {ujluj'i/ne){d^j^)o{dxxjz)o], 

= vi9yyjz)o + 2[{dyVy)o - {uj'Xuj'i/ne){dyiy)o]{dyB^)o 

+ 2Uj\dl[{dyVy)o{dyyi:,)o + (^y jy ) „ (5j/y V^ ) « 

- {ujluj'i/ne){dyiy)o{dyy],)o]. 

Substituting expressions f l23|) into these equations and using equations flTSjl . flT9|) and 
dxVx = —dyVy, we obtain 

- ri{dxx3z)o = 2{dyVy)o[{d^By)o - uj\dl{d^^jz)o] 

x[l+^{ujl-dl^/d^)l (24) 

-V{dyyJz)o = 2{dyVy)o[{dyB^)o + Uj\dl{dyyjz)o] 

x[l+7(u;!-42^M')], (25) 

where we introduce a useful dimensional parameter 

7 = ul{d,^yBz)o /ne{dyVy)o ■ (26) 

In the case of electron- ion plasma [ZrUe <^ rrii and w^ ~ w^ ~ 1), parameter 7 measures 
the relative strength of the Hall term (j x B)z/ne and the ideal MHD term (V x B)^ 
inside the electron layer. 

Taking the ratio of equations (j24l) and (!25l) . we obtain 

{dyB^)o = {dxBy)o{dyyiz)o/{dxx3z)o " 2W^rf^ ( 9yy j^ ) q 

^{B,,t6/L'){l + 2ujldll6'), (27) 

where we use the estimates {dxxjz)o ~ —jo/^^ and {dyyjz)o ~ —jo/L^, and equation ( l20l) . 
In equation fl2Tl) . the electric field E^ is balanced by the ideal MHD and Hall 
terms outside the electron layer, where the resistivity and electron inertia terms are 
insignificant. Therefore, 

Ez^ - V,By[l - {ujlujl/ne)3jV,] 

^{dyVy)o5Bext{l + ooll), (28) 

Ez ^ VyB^l - {uloj^^/ne)jy/Vy] 

~ {dyVy)o{dyB,)oL\l + c^^ ^) (29) 

at the points (x ~ 5, y = 0) and (x = 0,y ~ L) respectively. Here we use the estimates 
jx ~ {dxyBz)oS, jy ^ -{dxyBz)oL, 14 ~ -{dyVy)oS, Vy ^ {OyVy) oL , B^ ^ {dyBx)oL aud 
By ^ Bext, and equation fl26l) . The ratio of equations fl28|) and fl29|) gives 

(5y5,)o ^ 5e.i<5/L2 ^ BlJjoL\ (30) 

where we use equation (!20l) . Comparing this estimate with equation ( |27l) . we find 
(5 > cu+de ~ c^e- Therefore, using equation (l20i) . we obtain 

Jo < 5e.t/4 (31) 
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and Ez < rjBext/de [IS]- The estimate B^ ^ {dyBx)oL ~ Bext^/L for the value of 
the perpendicular magnetic field is in agreement with geometrical configuration of the 
magnetic field lines inside the electron layer of thickness S and length L. 
Combining equations ( 120|) . (!22|) and ( 128|) . we obtain 

vfo^{dyVy)oBUl + u'_^). (32) 

This equation describes conversion of the magnetic energy into Ohmic heat inside the 
electron layer with rate ^ |(c?xOo| =_l(f^xKr)o - ^l{dxjx)o/ne\ ^ {dyVy)o{l + 7) in the 
case of electron-ion plasma (wi ~ 1) I5II, and with rate ^ |(i9kK)o| = idyVy)o in the case 
of electron-positron plasma {u^ =0). 

Next, we use the 2;-component of Faraday's law, dxEy — dyEx = —dtBz = 0, where 
the time derivative is set to zero because we assume that the reconnection is quasi- 
stationary. We substitute Ex and Ey into this equation from Ohm's law ( !T7|) and, after 
tedious but straightforward derivations, we obtain 

Vidxjy - dyjx) + {uj'Xuj'i/ne){BxdxJz + Bydyj^) 

+ VxdxB, + VydyB, - BxdxV, - BydyV, 

~r ^j^O'^yVxKyxxJy (^xyjx) ~r VyyOxyJy Clyyjx) 

~r jxKyxx 'y (^xy 'xj T ]y\Oxy Vy ^yy ' xj 

- {uj'\uj^_/ne)jx{dxxiy - dxyjx) 

- iujlujl/ne)jy{dxyjy - dyyjx)] = 0. 

Taking the dxy derivative of this equation at the central point O and using equations (TTOll 
and f l23|) . we obtain 

0= -r] [{dxyxxBz)o + {dxyyyBz)o] + (^^ " dl^ / d^) 

X {uj1/ne)[{dxBy)o{dyyjz)o + {dyBx)o{dxx3z)o] 
^ rine{dyVy)ol/uj\5'^ 

- {ul - dl^/dl){ujl/ne)[jl/L^ + (5,i?J„jV5% (33) 

To derive the final expression, we use equation (!26|) and the estimates {dxyxxBz)o ~ 

-{dxyBz)o/S^ > idxyyyBz)o, {Oxxj z) o ~ -jo/ 5"^, idyyjz)o ~ -jo/L^, {,dxBy)o ^ jo- Uslug 

equations (IT5|) . f lTSj) . fl20|) and f l30|) . we rewrite equation ( 133|) as 

a;! - 427M' ^ r/L2(9,\/,),7K<Vi. (34) 

Note that equations ( 133|) and f l34l) result in 

< 7 < w-f^i/t^e- (35) 

Equation ( iT6l) for the plasma (ion) acceleration along the reconnection layer in the 
y-direction gives 

piVV)Vy + uldl{}V)jy = -dyP + J, 5, - J, S, . (36) 

^ In the case of electron-ion plasma, in the upstream region outside the electron layer the magnetic 
field lines are frozen into the electron fluid and inflow with the electron velocity u%. 
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Taking the y derivative of this equation at the central point O and using 
equations (IT5|) . f lT9|) and fl26|) . we obtain 

pidyVy^^il + dlf/dl) ^ BIJL^ + Jo{dyB,)o. (37) 

In the derivation of this equation we use the estimate {dyyP)o ~ {dyyBy/2)ext ~ 
—B'^^^/L'^, which reflects the fact that the pressure drop is approximately equal to 
the drop in the external magnetic field pressure. This estimate follows from the force 
balance condition for the slowly inflowing plasma across the layer, in analogy with the 
Sweet-Parker derivations □ [llj. Using equations ( IT8l) and ( l30l) . and neglecting factors 
of order unity, we rewrite equation ( 137|) as 

{dyVy)^^{VA/L){l+dl^yd^)-'/'. (38) 

Now we note that on the y-axis (x = 0) equation ( l36l) reduces to pVydyVy = 
—uj\dljydyiy " OyP + jz-Ba;- We Intcgrate this equation from the central point O 
to the downstream region outside of the ion layer, x = and y ~ Lext-, where 
ideal MHD applies and jy ~ 0. The plasma inertia term pVydyVy integrates to 
pVy /2 = {l/2){BextVy/VAY, the electron inertia term oj'^dljydyjy integrates to zero, 
the pressure term —dyP integrates to ~ B1^^, and the magnetic tension force term jzBx 
integrates to ~ Bl^^ Q As a result, we find that that the eventual plasma outflow 
velocity is approximately equal to the Alfven velocity, Vy ~ Va, in the downstream 
region outside of the ion layer (at y ~ L^xt)- 

In the end of this section, we derive an estimate for the ion layer half-thickness 
A. In these derivations we proceed as follows. Outside the electron layer the electron 
inertia and magnetic tension terms can be neglected in equation fl36|) . and we have 
p{yV)Vy ~ —dyP. Taking the y derivative of this equation at y = 0, we obtain 
p[Vx{dxyVy) + (dyVy)'^] ^ -{dyyP)o ~ B'^^JL'^. Hcrc the term Vx{dxyVy) is about of 
the same size as the term (dyVy) ^. Therefore, we find that {dyVy)ext ~ Va/L outside 
the electron layer (but inside the ion layer). Next, in the upstream region outside the 
ion layer ideal single-fluid MHD applies. Therefore, at x ~ A and y = equation fl2T|) 
reduces to E^ ^ -VxBy ^ -{dxVx)ext^Bext = {dyVy)ext^Bext ~ VAABext/L, where 
Ez is given by equation fl22|) . As a result, we obtain 

idyVy)ext ~ Va/L, a ^ r]joL/VAB,xt. (39) 

5. Solution for two-fluid reconnection 

To be specific, hereafter, unless otherwise stated, we will focus on two-fluid reconnection 
in electron-ion plasma and will assume Zmg <^ irti, de <^ di and w^ = w^ = 1. In this 

"*" For a proof, integrate equation (|16|) along the unclosed rectangular contour (x = 0, y = 0) — )■ (x = 
A, y = 0) — > (x = A, y = y) — > (a; = 0, y = y), then take the limit y — )■ and use the Taylor expansion 
in y for the physical quantities that enter equation ()16|) . For details refer to (llj . 

* Note that jz ~ jo for y < i and j^ ~ for y > L. Field Bx ~ {dyBx)oy ~ {B^xt / J oL^)v , see 
equation (|30l) . 
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2r,/j2 ^ ^r2/Q T/ \ ~ /j2t/2 y"^^) 



case equations (15^ and (IMIl reduce to 

r^Jo^ ^ idyVy)oBUl+^) 

l-dl^/d^ ^ VL\dyVy),^/dM 

We solve these equations and equations ( l20l) . ( I26l) . ( 130|) . ( l38l) and ( 139|) for unknown 
physical quantities jo, 5, A, L, 7, {dyVy)o, {dyBx)o and {dxyBz)o- We calculate the 
reconnection rate -E^ by using equation (!22l) . We neglect factors of order unity, and we 
treat the external field Bext and scale Lext as known parameters. Recall that parameter 
7, given by equation f l26|) . measures the relative strength of the Hall term and the 
ideal MHD term in the z-component of Ohm's law (in the case of electron-ion plasma). 
Depending on the value of parameter 7, we find the following reconnection regimes and 
the corresponding solutions for the reconnection rate. 

5.1. Slow Sweet-Parker reconnection 

When 7 < 1, both the Hall current and the electron inertia are negligible, the electrons 
and ions flow together, and the electron and ion layers have the same thickness and 
length. In this case, equations (138|1 and (HO!) become {dyVy)o ~ Va/L, rjj'^ ^ (^yK/)o-SLi 
and 1 ~ r]L'^{dyVy)o^/dfVA respectively. As a result, we obtain the Sweet-Parker 
solution [M[T51. 



1 « S = VAL^xt/v < LlJdl 

7 ^ VAd'i/7]Lext = Sdf/Ll^t, 

E^ ^ rj^/'Vl^'B^xt/Lli', = VAB,xt/S'l\ 

Jo ^ V^J^B^xth^'^Ll'xt = S^I^B,xtlL,xu 

6^A^ V'^'LIS/V}/' = L^xt/S'l^ > d,, (41) 

-^ ~ l^exty 

idyVy)o ^ idyVy)ext ~ VA/LexU 

{dyBx)o ^ v'^'B^xt/V^^'LlS = B,xtlLextS^I\ 
{dxyB,)o ^ VABextdi/r]Ll^^ = SB^xtdi/Ll^t, 

where the Lundquist number S ^ 1 is defined by equation (fT8l) . The condition 
S < Ll^^/df is obtained from 7 < 1. From this condition for S we find that Sweet- 
Parker reconnection takes place when di is less than the Sweet-Parker layer thickness, 
di < Lext/S^^'^, which is a result observed in numerical simulations |5l |6l [7]. Note that 
the quadrupole field is small in the Sweet-Parker reconnection case, B^ ~ {dxyBz)oL6 ~ 
{S^^'^di/ Lext)Bext ^ Bexti and the ion and electron outflow velocities are approximately 
equal to the Alfven velocity, Vy ^ {dyVy)oL ^ Va [SI E]- 

Now, let us for a moment consider the case of reconnection in electron-positron 
plasma. In this case d^ = di, ijj\ = 1/2, w^ = and equation fl35l) gives 7 = 0. This 
result represents an absence of the quadrupole field B^ [refer to equation fl2^ ]. which is 
known from numerical simulations [HI [T71 [18]. Therefore, our model predicts the slow 
Sweet-Parker reconnection solution for reconnection in electron-positron plasmas, which 
is in disagreement with the results of kinetic numerical simulations [T6l [T71 [T8] . A likely 
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reason for this discrepancy is that our model neglects pressure tensor anisotropy, which 
plays an important role in reconnection in electron-positron plasma. 

5.2. Transitional Hall reconnection 

When 1^7^ di/d^,, the Hall current is important but the electron inertia is negligible. 
In this case, equations ( l38ll and ( l40l) become {dyVy)o ~ Va/L, rjjl ^ {dyVy)oBl^ff 
and 1 ~ rjL'^{dyVy)o^ /dfVl. As a result, we obtain the following solution: 1 < 
7 ~ dfYAhL = Sdf/LL^^t < di/d^, E, ^ (di/L)Ki5e^t, jo ~ diVAB^^t/vL = 
SdiBext/LLext, S ^ rjL/diVA = LL^xt/Sdi, A ^ di, {dyVy)o ^ {dyVy)ext ~ Va/L, 
[OyBxjo ~ vBext/diVAL = B(,xtL(,xt/odiL, {OxyBzjo ~ diVABext/vL = odiBf^xt/ L L^xt- 
These results are in agreement with earlier theoretical findings [121 113 EOl El] . 

Condition 1 < 7 ^ di/d^ gives Sdedi/Lext < -^ < Sdj/Lext for the electron layer 
length L. Unfortunately, in our model, the exact value of L cannot be estimated 
in the Hall reconnection regime. In theoretical studies [I2l [191 El] length L was 
essentially treated as a fixed parameter. Here, we take a different approach and make 
a conjecture that the Hall reconnection regime describes a transition from the slow 
Sweet-Parker reconnection to the fast collisionless reconnection (presented in the next 
section). Numerical simulations and laboratory experiments have demonstrated that 
this transition happens when the ion inertial length is approximately equal to the Sweet- 
Parker layer thickness, di ~ Lgxt/VS O El [71 [221 ES]. Therefore, our conjecture leads 
to the following solution for the Hall reconnection regime: 

S = VALext/v ~ LlJdl 

Lext ^ L > deLext/di, 

7 ^ Lext/L, 

E, ^ idi/L)VABext, 

Jo '^ ^ext^extl (^i^i ( AO\ 

5 ^ diL/Lext > 4, 
A^d,>5, 

{dyVy)o ~ {dyVy)ext ~ Va/L, 

[OyBxjo ~ Bextdi/ LLext, 



2 



[OxyBzjo ~ B(,xtL(,xt/diL 

It is noteworthy that, in the Hall reconnection regime, the typical value of the quadrupole 
field is comparable to the reconnecting field value, B-^ ~ {dxyBz)oL5 ~ -Bext- The typical 
value of the ion outflow velocity is equal to the Alfven velocity, Vy ~ {dyVy)oL ^ Va- 
To estimate the typical value of the electron outflow velocity, we use equations (jSj), 
(dSD, (HHD and (IllD, and find ul ^ Vy - {m,/Zep)]y = Vy - {diVA/ Bext)3y ~ Ki + 
{diVA/B,xt){dxyB,)oL ^ VAiLext/L) > Va. 

As the electron layer length L decreases from its maximal value L fs L^xt to its 
minimal value L ?a deLext/di, the transitional Hall reconnection solution (H2|) changes 
from the slow Sweet-Parker solution (1411) to the fast collisionless reconnection solution 
presented below [see equations fH3l) - f[53l) and table [1]. 
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5.3. Fast collisionless reconnection 

When di/de < 7 < dfjdl [compare to equation (!35l)]. the electron inertia and the Hall 
current are important inside the electron layer and the ion layer respectively. In this 
case, equations (138!) and (140|) become {dyVy)ol ~ diVA/deL, rjj'^ ^ {dyVy)oBl^{) and 
1 — d'l'j/dj ~ riL'^{dyVy)o^/djVj^. As a result, taking into consideration equation ( IHT]) . 
we obtain the following solution: 

Lext/de < 5 = VALext/V ^ LJ^t/dedi, (43) 

cii/4 < 7 < <V^e, (44) 

E^ ^ r]Bext/de = iLext/Sde)VABext 

^ {A/L)VAB,,t ~ {d,/L)VAB,,t, (45) 

Jo ^ B,.,t/de, (46) 

5 ^ 4, (47) 

A ^ c^i > <5, (48) 

L ^ VAdedi/r] = Sdedi/Lexu (49) 

(5,K;)o ~ ^7/^27 = VALext/Sdl^ < Va/L, (50) 

(5yK,)e:rt ~ v/dedi = VAL^xtl SdJi ^ Va/L, (51) 

idyB,)o ^ B,,tr^^/VldJl = B^^tLlJS^dJl (52) 

{dxyBz)o ~ BextV/VAdldi = BextLext/Sdldi. (53) 

Here the limits on the Lundquist number given in equation (143|) . Lext/de ^ S < 
^Ixt/dedi, are obtained from the conditions E^ <^ Va-Bco;* (slow quasi-stationary 
reconnection) and L < Lgxt (the electron layer length cannot exceed the ion layer 
length). Except for the definition of the reconnecting field Bext, equations fHS]) - fH7|) 
and (H9|) essentially coincide with the results obtained in [13] for a model of electron 
MHD (EMHD) reconnection. The collisionless reconnection rate, given by equation ( 145|) . 
is much faster than the Sweet-Parker rate E^ ^ VABext/ ^~S [see equations (HTl)]. 

Note that the value of 7 or, alternatively, the value of the ion acceleration rate 
{pyVy)o ~ v/dll at the point O cannot be determined exactly. This is because in 
the plasma momentum equation (!36|) . the magnetic tension and pressure forces are 
balanced by the electron inertia term <ig(jV)jy inside the electron layer. The ion 
inertia term p(VV)V^ can be of the same order or smaller, resulting in the upper 
limit {dyVy)o < Va/L. In other words, inside the electron layer the magnetic energy 
is converted into the kinetic energy of the electrons (and into Ohmic heat), while the 
ion kinetic energy can be considerably smaller. Therefore, the ion outflow velocity 
can be significantly less than Va in the downstream region outside the electron layer 
(at y ^ L). At the same time, the electron outflow velocity is much larger than 
Va and is approximately equal to the electron Alfven velocity, Uy ~ {mi/Zep)jy = 
{diVA/ Bext){dxyBz)oL K. diVA/de ^ VeA = Bext/ y/nm^ > Va- Howcver, further in the 
downstream region, a.t y > L, as the electrons gradually decelerate, their kinetic energy 
is converted into the ion kinetic energy. As a result, the eventual ion outflow velocity 
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becomes ~ V^i, as was estimated in the end of Section |H These results emphasize the 
critical role that electron inertia plays in the plasma momentum equation flT6l) . These 
results also agree with simulations [27], which found the ion outflow velocity to be 
significantly less than Va_ in the downstream region outside of the electron layer, and 
found acceleration of ions further downstream (in the decelerating electron outflow jets). 

Our theoretical results for collisionless reconnection are in good agreement with 
numerical simulations and/or laboratory experiments!^. Indeed, the estimates A f« dj for 
the ion layer thickness, 6 ^ de for the electron layer thickness, Bz ~ {dxyBz)oSL ^ B^xt 
for the quadrupole field, and Uy ~ VeA = Bext/y/nm^ for the electron outflow velocity 
agree with simulations |Sl El 13 1251 ESI EZl [2H] • The estimates A ^ di and B^ ^ B^xt 
also agree with experiment [6]. However, the experimentally measured thickness of 
the electron layer is about eight times larger than our theoretical model and numerical 
simulations predict [29l|30]. This discrepancy can be due to three-dimensional geometry 
effects and plasma instabilities that may play an important role in the experiment [6l |30] . 

Our results are also in a qualitative agreement with recent numerical findings of 
an inner electron dissipation layer and of electron outflow jets that extend into the ion 
layer [23 [2S1 [23 [2H]- We note that the estimated electron layer length L ^ VAdedi/rj 
is generally much larger than both the electron layer thickness 6 ^ d^ and the ion 
layer thickness A ^ di, which is consistent with numerical simulations [23 |26l [27] . 
However, if resistivity r] becomes anomalous and considerably enhanced over the Spitzer 
value, then L can theoretically become of order of di and the reconnection rate can 
become comparable to the Alfven rate VA-Bext, which is also observed in numerical 
simulations [22| [28]. 

Unfortunately, a detailed quantitative comparison of our theoretical results to the 
results of kinetic numerical simulations is not possible because these simulations do not 
explicitly specify constant resistivity rj. In addition, in the simulations the anisotropy 
of the electron pressure tensor anisotropy was found to play an important role inside 
the electron layer and in the electron outflow jets [23 [2H]- In contrast, in the present 
study we assume an isotropic pressure, and the electrons are coupled to the field lines 
everywhere outside the electron layer (including the jets). 

In our model, the electric field Ez is supported by the Hall term (j x B)z/ne in the 
downstream region L <y < Lext- Therefore, in the collisionless reconnection regime, our 
model predicts an existence of Hall-MHD Petschek shocks that are attached to the two 
ends of the electron layer and separate the two electron outflow jets and the surrounding 
plasma. Note that, for electron-ion plasma {Znif. <^ rui), the ideal MHD and Hall terms 
in Ohm's law ( [T2l) can be combined together as —V x B + (mj/Zep)j x B = —u^ x B, 
where u'' is the electron velocity given by equation ([5|). Therefore, all results for the 
Hall-MHD Petschek shocks can be obtained from the corresponding results derived 
for the standard MHD Petschek shocks by replacing the plasma velocity V with the 



tt Even though reconnection rate (|45p is proportional to resistivity, we still use the standard term 
"collisionless reconnection" because in the fast reconnection regime rj should be viewed as the effective 
resistivity, which is to be calculated from the kinetic theory. 
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Table 1. Solution for two-fluid reconnection 





slow Sweet-Parker 


Hall 


fast 
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1 « 5 < LlJd] 


LlxJd^ 


Lext/de ^S < LlJdedi 
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Sdi /I/ea;t 


Lext/^ 


d^/de < 7 < dj/dl 


E, 


VA5e.,/Sl/2 


{d,lL)VAB,,t 


{Lext/ Sde)VABext 
^ {djL)VABext 


Jo 


•S* Bext/Lext 


BextLext/diL 


Bext/de 


6 


Lext/S'/' > d. 


diL/Lext > 4 


de 


A 


L,xt/S"^ ^ 5 


d^>S 


di-^S 


L 


J-'ext 


Lext ^ L > df,Lf.xt/di 


Sdedi/Lext 


idyVy)o 


Va/L 


Va/L 


VALext/Sdl^ < Va/L 


{dyVy)ext 


Va/L 


Va/L 


VALext/Sded, ^ Va/L 


idyB,)o 


Bext/LextS 


Bextdi/ LLext 


BextL^^JS dedi 


{dxyB,)o 


SBextdi/L^^^ 


BextLext/diL 


BextLext/ bd^di 



electron velocity u''. In particular, the parallel components of the magnetic field and 
electron velocity jump across the Hall-MHD Petschek shocks, the velocity of the shocks 
is K. \ul\ ^ {mi/Zep)\i^\ ^ {diVA/Bext){dxyB^)o5 ^ VALext/Sde < Va, and the opening 
angle between the shocks is ~ B^/By ^ {dyBx)oL/Bext ~ Lext/Sdi <^ 1. Shocks were 
indeed observed in numerical simulations j31j . However, in these simulations a spatially 
localized anomalous resistivity was prescribed, resulting in a short layer length, while 
in our study resistivity r] is assumed to be constant. 



6. Discussion 

The solution for two-fluid reconnection is summarized in table [H This table includes 
solution formulas for three reconnection regimes: the slow Sweet-Parker reconnection 
regime, the transitional Hall reconnection regime, and the fast coUisionless reconnection 
regime. The reconnection rates for these three regimes are respectively shown by the 
solid, dotted and dashed lines in figure O 

It is well known that resistivity rj can be considerably enhanced by current-driven 
plasma instabilities [HI [71 [21]. Because the coUisionless reconnection rate Ez ~ rjBext/de 
is proportional to the resistivity [see equation fHSl) ]. this rate can increase significantly 
as well. As a result, we propose the following possible theoretical explanation for the 
two-stage reconnection behavior (fast and slow) that is frequently observed in cosmic 
and laboratory plasma systems undergoing reconnection processes. 

During the first stage, such a system is in the very slow Sweet-Parker reconnection 
regime, during which magnetic energy is slowly built up and stored in the system. The 
magnetic energy and electric currents build up, the field strength increases and the 
resistivity decreases [32]. As a result, the Lundquist number S increases and the system 
moves to the right along the solid line in figure [21 
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Ez \S-' 


di/de 




1 


»A 

1 »► 



Lext/S = dp S 

Figure 2. Schematic plot of tlie reconnection rate E^ versus tire Lundquist number 
S in the slow Sweet-Parker (solid line), transitional Hall (dotted line), and fast 
collisionless (dashed line) reconnection regimes. 



When the Lundquist number S becomes comparable to L^^^/dj and the thickness 
of the current layer L^.^t/S^^'^ becomes comparable to di, the system reaches point A 
in figure O Next the system goes into the transitional Hall reconnection regime and 
quickly moves up along the vertical dotted line in figure [2j During this transition, 
the length of the electron layer shrinks from ^ L^^t to ~ {de/di)Lext, the electron 
layer thickness decreases from ^ di to ^ de, and both the electric current and the 
reconnection rate increase by a factor ^ di/de ^ 1. The system ends up in the fast 
collisionless reconnection regime at point B in figure El 

Because of the considerable increase in the electric current during the Hall 
reconnection transition from point A to point B, plasma instabilities develop, and, 
consequently, resistivity rj becomes anomalous and rises in value. As a result, the 
reconnection rate E^ ~ rjBext/de increases, the Lundquist number S = VALext/v and 
electron layer length L ^ VAdedi/r] decrease, and the system moves from point B to the 
left along the dashed line in figure [2j The system enters the second stage characterized by 
a rapid release of the accumulated magnetic energy. Even though our theoretical model 
is stationary, assumes constant resistivity and cannot describe this stage in detail, the 
physical mechanism of slow and fast reconnection outlined above is self-consistent and 
may take place in nature. 
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